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Abstract 



. We introduce a method to construct wave packets with complete classical and quantum corre- 

spondence in one-dimensional non-relativistic quantum mechanics. First, we consider two similar 
oscillators with equal total energy. In classical domain, we can easily solve this model and obtain the 
trajectories in the space of variables. This picture in the quantum level is equivalent with a hyper- 
bolic partial differential equation which gives us a freedom for choosing the initial wave function and 
its initial slope. By taking advantage of this freedom, we propose a method to choose an appropriate 
initial condition which is independent from the form of the oscillators. We then construct the wave 
packets for some cases and show that these wave packets closely follow the whole classical trajecto- 
ries and peak on them. Moreover, we use de-Broglie Bohm interpretation of quantum mechanics to 
quantify this correspondence and show that the resulting Bohmian trajectories are also in a complete 
O^l , agreement with their classical counterparts. 
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1 Introduction 

The issue of classical-quantum correspondence has been extensively investigated in the literature pQ. 
Moreover, the question of construction and interpretation of wave packets in quantum mechanics and 
its connection with classical mechanics has been attracting much attention. In quantum physics, one is 
generally concerned with the construction of wave packets by the superposition of the energy eigenstates 
which would peak around the classical trajectories. These efforts have been started by Schrodinger [5] 
and followed by others who were interested in finding quantum mechanical states which provide a close 
connection between classical and quantum formulations of a given physical system in the context of the 
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coherent states [3|. These states can also be generated using algebraic methods [3], supersymmetric 
quantum mechanics [5] and its application to different physical situations [5J. 

In this paper, we pursue a different approach to construct wave packets with complete classical and 
quantum correspondence. First, we consider similar Schrodinger equations of two oscillators with u and 
v variables. Classically, we can easily solve this model and find the behavior of the variables versus 
time (u(t),v(t)). If these two oscillators have the same total energy, the solutions will be equivalent 
up to a temporal phase factor. Imposing particular initial conditions, we will obtain the trajectories 
in the configuration space (u, v). Inversely, using the form of the trajectories, again we can find the 
temporal behavior of the variables. For instance, in general, the trajectory for the case of a simple 
harmonic oscillator (SHO) is an ellipse. Therefore, we can obtain the time behavior of the variables 
as u(t) — Acos(ujt) and v(t) — Asm(ujt + </>g), where 4>q = corresponds to a circle. Thus, we can 
parameterize the behavior of one variable (e.g. u) versus another variable (e.g. v) where they both 
satisfy the same equation of motion. 

Our motivation for this model is related to its interesting quantum mechanical properties. In the 
quantum version, this model allows us to define a wave function which is a function of u and v unlike the 
usual case where the wave function is a function of u or v and time. As we shall see, the existence of the 
two mentioned similar oscillators results in a hyperbolic partial differential equation (PDE). Therefore, 
we are free to choose the initial wave function and the initial slope of the wave function. This freedom 
for choosing the initial conditions (or the expansion coefficients) allows us to obtain a specific form of 
the wave packets with complete classical and quantum correspondence. In fact, the initial wave function 
and its derivative correspond to classical initial position and initial momentum, respectively. To fix the 
expansion coefficients, we will apply the prescription used in Refs. 7,8,9 for the case of hyperbolic 
PDEs which also appear in the context of quantum cosmology. 

The purpose of this paper is to construct wave packets such that the classical and quantum corre- 
spondence is manifest, which means that the wave packet should be centered around the classical path 
and the crest of the wave packet should follow the classical path as closely as possible. Usually, the 
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temporal behavior of the variables is studied in literature where its underlying equation is parabolic. 
But here, we are encountered with the hyperbolic PDEs which correspond to the classical picture of 
parametric trajectories. Although both descriptions are equivalent, the latter contains some interesting 
properties in quantum domain which allows us to obtain the complete classical-quantum correspondence 
for any oscillator using a unique prescription. 

To be more precise, we can also quantify this correspondence using de-Broglie Bohm interpretation of 
quantum mechanics [lOj . This approach gives us the Bohmian trajectories which are guided by the wave 
function. These trajectories are governed by classical plus quantum potential. Therefore, the coincidence 
between classical and quantum trajectories accrues in the limit of vanishing quantum potential. In fact, 
the usage of the appropriate initial conditions results in the suppression of the quantum potential [5] . 

The paper is organized as follows: In Sec. [2l we present the model of two similar oscillators which 
results in a hyperbolic PDE at the quantum level. We then construct wave packets by choosing the 
appropriate expansion coefficients which leads to a good correspondence between classical and quantum 
solutions. In Sec. [H we use the casual interpretation of quantum mechanics to quantify this correspon- 
dence and obtain the Bohmian trajectories for different cases. In Sec. IH we state our conclusions. 



In non-relativistic quantum mechanics we can write the dimensionless time independent Schrodinger 
equation for variable u as 



where V is the potential term which here is supposed to be an even function of its variable and E is the 
total energy. Also we can rewrite the above equation for another variable v 



2 The model 



V(u) 
du 2 



+ V{u)*/j(u) = Ei()(u), 



(1) 



dv 2 



+ V(v)ip(v) =Eip(v). 



(2) 
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Classically, the above equations show two similar oscillatory motions for u and v variables with equal 
total energy and the following equations of motion 

dV(u) .. dV(v) 

By assuming specific classical initial conditions, we can obtain a particular trajectory in the configuration 
space. Here, we impose the following initial conditions: 

w(0) = 0, tt(Q)=«o, «(0)=tib, u(0) = 0, (4) 

which result in circular motion for the case of simple harmonic oscillators for u and v variables (Fig. [2]). 
Moreover, the classical behavior of u versus time in the presence of some sample potentials is also 
depicted in the right part of Figs. [3][6] as dashed lines. Note that, the above initial conditions result in 
classical solutions which have the following property 

v(t) = u(t - T/4), (5) 

where T is the classical period of motion. 

Now, we define a new wave function ^{u, v) = ip(u)ip(v) where ip(u) an d V'( w ) satisfy Eqs. ([I} and 
([2|), respectively. Therefore, *f?(u,v) satisfies the following partial differential equations 



n du (6) 



t V(v) *(it, v)=E v). 



Subtracting these equations leads to 



where is a hyperbolic differential equation. This equation can be solved using the separation of variables 
and the general wave packet which satisfies this equation can be written as 

*(«,«)= ^ A n il} n (u)ip n (v) +i ^ B n i/j n (u)^ n (v). (8) 

n=even n=odd 

Since the potential is an even function of its variable, the eigenfunctions are separated into even and odd 
categories. Moreover, we have intentionally separated the odd and even terms for further usages. To 



find the exact form of the wave function we need to specify the expansion coefficients. These coefficients 
will be determined from the initial form of the wave function which is evaluated along u— axis (v = 0) 
in consistency with the classical initial conditions ((4]). 

Now, let us consider the initial behavior of the wave function. The wave function and the first 
derivative of the wave function at v = take the form 



*(«,0)= A n t(> n (u)ip n (0), (9) 



d^(u, v) 



= i 

v =° n=odd 



J2 Bn^«(«X(0)- (1°) 



dv 

Therefore, the coefficients A n determine the initial wave function and the coefficients B n determine 
the initial derivative of the wave function. From a mathematical point of view, since the underling 
differential equation Q is second order, A n s and B n s are arbitrary and independent variables. On the 
other hand, if we are interested to construct wave packets which simulate the classical behavior with 
known classical positions and momentums, all of these coefficients will not be independent. It is obvious 
that the presence of odd terms dose not have any effect on the form of the initial wave function but 
they are responsible for the slope of the wave function at v = 0, and vice versa for the even terms. Now, 
consider the behavior of the initial wave function. Near v — the differential equation takes the 
form 

{-^2+^2+ V(u) - F(0)} ^(u,v) = 0, (11) 

which has the solution as 

ip(u,v) =ip(u)x(v), (12) 

where ip( u ) an d x( v ) satisfy 

d 2 Xn(v) , , ( . 

—fop. yE n x n (v) = 0, (13) 

~ Z 2 +(V(u)-V(0))^ n (u) = E n ip n (u), (14) 

and E^s are separation constants. These equations are Schrodingcr equations with En's as their energy 



levels. Equation (fT5|) is exactly solvable with plane wave solutions 

Xn(v) — a n cos C\/E^ + iP n sin (^/Kvj, (15) 

where a n and (3 n are arbitrary complex numbers. We can find the eigenfunctions and the eigenvalues of 
equation (|14|) using numerical techniques like Spectral Method [11) with an acceptable accuracy. Now, 
the general solution to Eq. (fTTjl can be written as 

ip(u,v) = ^ A* n cos(\f~E^v)il) n {u) + i ^ B n s\n(\/E^v)ip n {u). (16) 

n=even n=odd 

As stated before, this solution is valid only for small v. Therefore, we can write the initial conditions as 
follows 

lf>(tt,0) = ^4*>n(u), (17) 

even 

V-'KO) = iJ2K^Mu), (18) 

odd 

where prime denotes the derivative with respect to v. Obviously, a complete description of the problem 
would include the specification of both of these quantities. However, since we are interested to construct 
wave packet with classical properties, we need to assume a specific relationship between these coefficients. 
The prescription is that the coefficients have the same functional form [H [9] i.e. 

A* n = C{n) for n even, (19) 
B* = C{n) for n odd, (20) 

where C(n) is a function of n. In terms of A n s and B n s we have 

C(n) for n even, (21) 



B " = TH^ C ( n ) for n odd. (22) 
First, let us apply this method to the simplest (bounded) case which is an infinite square well 

3 -#<g<f, 

Vlq)={ (23) 
x) otherwise, 



B n = \j — C(n) for odd terms, (26) 



with well-known orthonormal even and odd eigenfunctions (here q stands for u or v) 

/f cos( !i££) n= 1,3,5,..., 

rj.q) = { (24) 
^sin(^) n = 2,4,6,.... 

2 2 

Since the energy spectrum for this model is E n — n ^ , using the exact form of the eigenstates, we can 
find the desired expansion coefficients (|2"Tj) 

A„ = y^C(n) for even terms, (25) 

2 

which result in the following wave packet 

*( u ,i;) = W-^ C(n)cos( — )cos( — )+t C(n)sm( — )sm( — ) L (27) 

ln=l,3,5,... n=2,4,6,... J 

To completely determine the wave packet, we need to specify the functional form of C{n). On the other 
hand, we can choose C(n) by imposing an appropriate initial condition (|19p which we choose as two 
Gaussian peaked at u = ±d 

#(u,0) = e - Q ("- d ) 2 + e -*{u+d)\ (28) 
This choice of initial condition is equal to the following form of the expansion coefficient 



-ze 4«l^ 

C(n) = 



v / 2aL/n 



_ , ri7r + 2iaL (d — L) . 

-ern — 

V 2^^ ; 



, . mr — 2ictL (d — L) . , nit — 2iaL (d-\- L) . 

- e L l erfi( — — } - crfi( — vsz — } 

„ , n7r + 2iaL (d + L) . 

crfil = I 

V 2v^L ' 



(29) 



where erfi(<7) is the imaginary error function erfi(q) = — ierf(iq). In Fig. [1] we have shown the resultant 
wave packets for different values of d and a. Classically, the particle is free inside the well and can begin 
its one-dimensional motion with positive or negative initial velocity and from any position between 
±L/2. As it can seen from the figure, the behavior of the wave packets are in complete agreement 
with the classical scenario and they strongly peak on the classical trajectories. In fact, the value of d 
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Figure 1: Infinite square well: the upper view of the square of the wave packet \il)(u,v)\ 2 for d = 3.5, 
a = 10 (left), d — 2.5, a — 10 (middle), d = 1.5, a — 5 (right) and L — 10 for all cases. 

determines the initial classical position. The presence of two rectangles with opposite direction shows 
the two possible direction of initial motion at u = ±d. Since the wave packets follow some straight 
lines, its classical picture consists of a free particle with constant velocity. Moreover, the height of the 
crest of the wave packet is constant along the classical trajectories which shows that the probability of 
finding the particle is constant along the classical path. On the other hand, since the classical velocity 
of the particle for this case is a constant of motion, the classical probability of finding the particle is also 
constant along its trajectory. Although this method is generally applicable for one-dimensional systems, 
we can also obtain a class of two-dimensional classical and quantum correspondences for this especial 
case. As the figure shows, these wave packets correspond to the cases of equal absolute velocities in u 
and v directions (|u u | — \ v v\) with arbitrary initial position. 

Now, we study other cases with polynomial or exponential form of the potential. First, we should 
choose an appropriate functional form of C{n). This means that, we need to specify C(n) in such a way 
that the initial wave function has a desired classical description (see Figs. EE]). We will see that this 
choice of coefficients leads to complete classical and quantum correspondence in all the studied cases. In 
the left part of Figs. EHU we have depicted the wave packets for some sample potentials. As it can be seen 
from the figures, the wave packets closely follow the classical paths and peak on them. Before discussing 
the solutions, we outline the causal interpretation of quantum mechanics and use it to interpret the 
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results in the next section. 

3 The causal interpretation 

To make the connection between the classical and quantum results more concrete, we use the causal 
interpretation of quantum mechanics [5J IIP) . In this interpretation the wave function can be written 
as 

= Re iS , (30) 
where R = R{u, v) and S = S(u, v) are real functions and satisfy the following equations 



du 2 dv 2 \du ) \ dv 

R d 2 S R d 2 S 2 dRdS_ ^dRdS^ _ 
du 2 dv 2 du du dv dv 

To write R and S, it is more appropriate to separate the real and imaginary parts of the wave function 

^(u,v) — x(u,v) + iy(u,v), (33) 
where x,y are real functions of u and v. Using equation (|30p we have 



R = ^x 2 +y 2 , (34) 

S = arctan( — ). (35) 

x 

On the other hand, the Bohmian trajectories are governed by 

Pu= ~du' ^ ' 

Vv = ^, (37) 

where the momenta correspond to the classical related Lagrangian L{q) = q 2 - V(q). Therefore, the 
Bohmian equations of motion take the form 

. i i d /r 



, , ( y ~) , (38) 
2 i + (1) du \xJ 

2 1 + (i) 2 ^(l)' (39) 



These differential equations can be solved numerically to find the time evolution of u and v. In the right 
part of Figs. we have shown the Bohmian trajectories as solid lines for five different forms of the 
potentials 

f V 1 (q) = q 2 , 

V 3 (q) = -3 9 2 + q\ (40) 
V 4 (q) = exp(q 2 /8) - 1, 
V 5 {q) = cosh(q) - 1. 

The initial wave functions are chosen to be localized on classical initial positions and also the initial slope 
is determined by equations (|17[) and (|19[) . In fact, we are free to choose any arbitrary but appropriate 
initial wave functions which correspond to the classical scenario. For cases studied here, we choose two 
different forms of coefficients (Figs. [2JS]) 

C (n) = _^ e -C 2 /4 c{n) = J^ e -e/\ (41) 

These coefficients are chosen so that the initial wave function contains two bumps along u— axis which 
correspond to initial and final classical positions. In fact, these choices are not special and there are 
many other equivalent alternatives. Since the slope of the wave function is related to the odd terms, 
the expectation value of the momentum operator is positive and negative respectively for initial and 
final states of motion in agreement with the classical scenario. Note that, since these two bumps are 
related to two different physical situations, the integrations are along positive and negative parts of 
the u — axis, respectively. Figure [5] shows the resulting wave packet for the simple harmonic oscillator 
and the Bohmian and classical trajectories. In the right part of figures [3][6l we have also shown the 
Bohmian and classical behaviors of variable u versus time (u(t)). The complete correspondence between 
these results are again manifest. It is worth to mention that for the cases where u(0) ^ 0, we need to 
choose complex C(n) [12] . Beside the stated advantages, this method has also some limitations. As 
we have shown by various examples, this method works well for symmetric bounded potentials. This 
class of potentials results in oscillatory classical motions which can be followed by wave packets with 
a finite set of eigenfunctions. Therefore, this method, in the present fashion, is not applicable even for 
simple unbounded situations like free particle, potential barrier, delta function potential and etc. The 
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generalization of this method for the case of free particles is the subject of our future work. 

4 Conclusions 

In this work, we have studied the classical-quantum correspondence in the context of non-relativistic 
quantum mechanics. First, we considered two similar oscillators for u and v variables and obtained 
the classical trajectories. Since these two variables satisfy the same equation of motion, we obtained 
the trajectory of one variable versus another one by choosing a particular set of initial condition. In 
the quantum version, this scenario leads to a hyperbolic partial differential equation where its solution 
contains infinite unknown coefficients. Half number of these coefficients arc related to the initial form 
of the wave function and others correspond to the initial slope of the wave function. Upon using 
a specific relation between these coefficients and choosing appropriate initial wave functions we can 
construct desired wave packets. The crests of these wave packets closely follow the classical paths from 
initial position to the final position. We have quantified this correspondence using de-Broglie Bohm 
interpretation of quantum mechanics. We applied this method to various cases and showed that the 
resulting wave packets completely simulate their classical counterpart's behavior. In particular, the 
Bohmian trajectories coincided well with classical trajectories for the all cases. 
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Figure 4: Double- Well potential: Left, the square of the wave packet \ip(u,v)\ 2 for C(n) 
and £ = 3. Right, the classical (dashed line) and Bohmian (solid line) trajectories. 
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Figure 5: V±{q) — exp(q 2 /8) — 1: Left, the square of the wave packet \ip(u, v)\ 2 for C(n) 
and £ = 3. Right, the classical (dashed line) and Bohmian (solid line) trajectories. 
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Figure 6: Vs(g) = cosh(g) — 1: Left, the square of the wave packet \ip(u,v)\ 2 for C(n) = —^L==e ^Z 4 
and C = 4. Right, the classical (dashed line) and Bohmian (solid line) trajectories. 
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